Relaxation of flying spin qubits in quantum wires by hyperfine interaction 
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We consider the relaxation of a spin qubit in a quantum dot propagating as a whole in a one- 
dimensional semiconductor with hyperfine coupling. We show that this motion leads to qualitatively 
new features in this process compared to static quantum dots. For a fast straightforward motion, 
the initial spin density decreases to zero with the relaxation rate independent of the spatial spread 
of the electron wave function and inversely proportional to the electron speed. However, for the 
oscillatory motion, the qubit acquires memory, and the dephasing becomes Gaussian rather than 
exponential. After some time, one third of the initial spin polarization is restored, as it happens 
for static qubits. This revival can occur either through periodic peaks or through a monotonous 
increase in the polarization, after a minimum, until a plateau has been reached. Our results can be 
useful for the understanding of the spin dynamics and decoherence in quantum wires. 

PACS numbers: 72.25.Rb,73.21.Hb,31.30.Gs,03.67.-a 



Semiconductor quantum wires play an important role 
in fundamental and applied spintronics. They allow one 
to study a variety of interesting phenomena related to 
spin transport [J- 11| and open the way to realizations of 
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exotic states such as Major ana fermions 

One of the most interesting aspects of spintronics is the 
hyperfine coupling of spins of carriers to the nuclear mag- 
netic moments of the host lattice 14 - 19| , strongly depen- 
dent on the spatial extension of the carriers wave func- 
tions. In semiconductor nanostructures, it is enhanced 
compared to the bulk due to confined wave functions in 
quantum dots [2ol422l and due to localization by disorder 
in quantum wires. |23| The spin dephasing times observed 
in quantum dots [2l| agree very well with what can be 
expected from the theory. [20| 

The observed relaxation times suggest that a single 
electron spin could become the physical realization of a 
quantum bit (qubit), crucial for quantum information 
applications. Unlike the dots, where electrons are well 
localized, extended nanowires are natural channels for 
traveling qubits, hereafter referred to as "flying" qubits 
3, [HI- Spin decoherence is the major concern for using 
qubits in information processing. At low temperatures, 
it is mainly caused by spin-orbit and hyperfine interac- 
tions. In undoped wires, the mechanisms for spin-orbit 
relaxation can be canceled out by tuning their cross sec- 
tion geometries for the wires extended along the [001] 
direction (Dresselhaus term) and also by keeping their 
environments symmetric (Rashba term). Hyperfine cou- 
plings hence can become the only source of spin depolar- 
ization. Their influence on the spin dynamics of qubits 
propagating through nanowires should be deeply under- 
stood for spin transport devices. Recently, Huang and Hu 
(26| have proposed a formulation of the spin-relaxation 
problem well suited for quantum information purposes. 
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FIG. 1: (Color online) Schematic plot of a semiconductor 
quantum wire. Random arrows show the magnetic moments 
of nuclei. Dashed lines correspond to density distributions 
with given expectation values of the z coordinate at differ- 
ent times: For n, of the order of T£, the matching of tails 
yields some correlation; for T2, considerably larger than r^, 
correlation is negligible. 



This formulation comes from studying the decoherence 
of an electron spin that is being transported in a random 
electric potential, as enclosed in a quantum dot with the 
time-independent shape of the wave function. The ability 
of high-fidelity electron transfer between quantum dots 
by surface acoustic waves was demonstrated recently. [27l — 
[30[ Here we apply a similar approach to study the effects 
of hyperfine interaction on spins in moving quantum dots, 
and show that the dynamics brings about new time scales 
and qualitative features compared to the static dots. 

We consider nanowires with L x L cross section, see 
Fig. [TJ and we take GaAs as a representative spintron- 
ics material. Nuclear magnetic moments are considered 
frozen on the time scale of the electron precession [20| , 
and are represented as random 1 = 3/2 spins with the 
expectation values of components Ii, where i = x,y,z. 
Under the assumption of decoupled spin and orbital mo- 
tion, valid in the absence of spin-orbit coupling and for 
weak hyperfine interaction, we can factorize the total 
wave function as ip{r,t)\^>). Here the orbital state is 
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^(r,£) and the spinor state is \<p) . The resulting Hamil- 
tonian for the spin degree of freedom depends on the 
orbital wave function as follows: 
Av ^ 



hf : 



— ^ |#„,t)| 2 (/W^ + I^cjy + I^a z ), (1) 



where the summation goes over all nuclei with spin com- 
ponents ij n ' at sites r n , A= 45 /ieV is the typical hyper- 
fine coupling constant for GaAs [io| , v is the volume per 
single nucleus (one eighth of the unit cell) , and oi are the 
Pauli matrices. 

Hamiltonian (pQ) can be rewritten in the form 

H M =\{n.v), (2) 

where the precession due to hyperfine interaction ft is an 
integral that depends on the electron position and can be 
expressed as 

«i = X J \^,t)\ 2 Pi (r)dV. (3) 

Here the density of the i-th component of the nuclear 
magnetic moment Pi(r,t) satisfies the white- noise distri- 
bution ((pi(r)pi(r f ))} = poS (r — r') , where ((. . .)) stands 
for the ensemble average and po = (lf)/v; here v~ x is 
the concentration of nuclei. To describe a moving elec- 
tron, we need the characteristic time-dependent correla- 
tor which can be obtained, assuming Gaussian random 
III, [II, 



fluctuations 
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«a(t)fii(t + r)» = 
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\^{v U t)\^{Y 2 ,t + T)\ 2 (( Pi {v 1 ) Pi {v 2 )))dV 1 dV 2 

= P0 J J \^(r,t)\ 2 \^(r,l 

Although, due to the isotropy of the random distribution 
of nuclear spins, these correlators are i— independent, 
for definiteness we consider the z— component of the 
spin below. With the course of time, the overlap of 
the electron wave function with itself decreases and the 
correlation gradually vanishes, as can be realized from 
FigHJ The i/;(r n ,t) states characterizing the propagating 
qubits are defined as i/j(r n ,t) = T](z n ,t)ip(x n ,y n ), where 
(f(xn,yn) = s'm(7rx n /L) sm(7ry n /L) describes an 

electron in a single- mode box, and rj(z n ,t) is assumed to 
be a Gaussian centered at the time-dependent position 
(z(t))j corresponding to the ground state of the moving 
parabolic potential: 

V(zn) = (-!^) V4 exp [-(*„ - (z(t)}) 2 /^ 2 } ■ (5) 

For this choice of the wave function, the correlator in 
Eq.((4j) is expressed as 

((fii(t)fii(f + r)» = exp [- (Az(t)) 2 /4£ 2 ] , (6) 
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FIG. 2: (Color online) Spin dynamics of static qubits for dif- 
ferent width £ as solved by 4 th — order Runge-Kutta method 
and averaged over 4096 realizations of disorder in the nuclear 
magnetic moments with the statistical errors of the order of 
1/V4096 = 1.6%. Here and below we assume the wire width 
L = 8 unit cells. For this wire width and £ = 10 nm, we ob- 
tain 1/ ((Hf )) ~ 1.5 ns, corresponding well to the relaxation 
time in the Figure. Since the GaAs lattice constant is 0.565 
nm, the nanowire cross-section is 4.5 x 4.5 ~ 20 nm 2 . 



where ((ft?)) = 3 (A/hf 1(1 + l)is/8^L 2 £ and Az(r) = 
(z(t + r)) — (z(t)). This correlator reaches the maximum 
when the electron is back to the position occupied at time 
t. For motion with constant speed, Az(r) = vr and the 
corresponding effective correlation time can be defined as 

Now we analyze the evolution of polarization of static 
and propagating wave packets. The spin dynamics is 
calculated for ensembles of random wires and averaged 
afterwards. At the initial time, qubits are fully polarized 
in the nonequilibrium spin state 1 0(0) ) = |1). Later on, 
for t > 0, | </>(£)) evolves following the time dependent 
Schrodinger equation, 
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d\ 



dt 



H M {t) 



(7) 



where the time dependence of Hht(t) is due to the dis- 
placement of the qubit. The resulting spinor \4>{t)) is 
related to the initial spin state as 



\4>(t)) = Texp 



H M {t)dt 



(8) 



where T stands for the time ordering. 

We consider three different realizations of spin dephas- 
ing: for static qubits (to have the known system to com- 
pare with), and for straightforward and oscillatory mo- 
tion. When packets are fixed to initial positions, the 
calculated average polarization decays to a minimum of 
about 0.1 of the initial value and monotonically increases 
up to a constant plateau where ((s z (t))) = ((s z (0)))/3. 
We consider packets with widths of £ = 5, 10, and 20 
nm; their spin evolutions, plotted in Fig. [2j look exactly 
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FIG. 3: (Color online) Average polarization for qubits in 
straightforward motion for different speeds v and widths £. 



as expected from the dynamics induced by hyperfine cou- 
pling in quantum dots [lo|, [HI . This steady polarization 
arises due to the precessional motion of the qubit spin 
around the time-independent effective field of the nuclear 
magnetic moments [20j. The spin evolution presented in 
this Figure is universal in the sense that it can be fully 
described by two timescales originated from the same 
random position-dependent spin precession: that for the 
initial polarization to decay down to 1/e, and that for 
reaching the 1/3 plateau. The scale for the initial de- 
cay rate of static electrons is then r st = \J ((Of )), and 
the entire time dependence in Fig. [2] can be described by 
using this single parameter [2o[. The dependence of the 
spin dynamics on the qubit spatial width is related to the 
number of interacting nuclei inside the electron cloud: A 
smaller amount of nuclei yields stronger fluctuations in 
the field resulting in a faster decay. 

In the next, straightforward motion regime, packets 
propagate along the wire with constant velocities v = 10 
and 20 nm/ns, taken here as examples. A behavior con- 
siderably different from the static one emerges when r^T st 
is much less than 1, that is, when the electron moves to 
a strongly different realization of the hyperfine field be- 
fore sufficiently changing its spin; this behavior is shown 
in Fig. [3l To gain insight into the influence of width 
and speed, we have taken into account, again, three dis- 
tinct packet sizes, £ = 5, 10, and 20 nm. When v = 10 
nm/ns all polarizations decrease below 0.1 after 10 ns and 
monotonically tend to zero afterwards. In addition, the 
smaller the qubit width, the faster the dephasing. Sim- 
ilar calculations for v = 20 nm/ns show a considerably 
less efficient decoherence. This behavior can be detailed 
as follows. The average spin evolution can be described 
with 13211 
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For the exponentially decaying part of the dynamics, this 
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FIG. 4: (Color online) Average polarization for qubits flying 
back and forth along the wires for various speeds and periods, 
£ = 10 nm. Note the sharp revivals in the regime when T = 30 
ns and v = 10 nm/ns. 



expression has the form 

«**(*)» 



exp(-y/TTT v t), 



«**(<>)» 

where the spin relaxation rate is determined by 

r„ = «^ 2 )) T? . 



(10) 
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Therefore, the fast motion establishes a new and smaller 
relaxation rate, T v = T st • (r st r^), such as T v <C r st . 
Here T v does not depend on the width of the electron 
wave packet but only on its speed, as we explain be- 
low using a scaling argument. The electron spin preces- 
sion rate, due to the frozen nuclear spins, is of the order 
of O ~ (A/%) /y/Nj where the number of nuclear mag- 
netic moments per single electron is N ~ L 2 ^/v. From 
Eq. ([TT]) . the resulting T v is proportional to (A/%) 2 v/L 2 v 
and, henceforth, also £— independent. This argument can 
easily be applied to other spatial dimensions D in the fol- 
lowing way. The square of the spin precession rate being 
of the order of ((Of)) - N' 1 ~ £~ D , while r - £/v, 
is D— independent. As a consequence, T v ~ £ 1_jD /^, 
which tends to zero for bulk crystals (D = 3) and two- 
dimensional structures (D = 2) when the size of the wave 
packet is increased. However, it remains constant in the 
wires (D = 1), making the results only weakly dependent 
on the details of the electron wave function, provided that 
it broadens relatively slow, on the timescale of spin pre- 
cession in the hyperfine field, such that the fast motion 
condition r^T st <C 1 holds. We notice that, for the elec- 
trons transferred by the surface acoustic waves 27-30| 
with the speed close to 3 x 10 3 nm/ns, the condition of 
fast motion is well satisfied for £ up to 500 nm. As a re- 
sult, the spin states are well protected against relaxation. 

A significantly different behavior is obtained when the 
packet moves back and forth, in a regime which can 
be experimentally achieved by the surface acoustic wave 



technique. [29j We analyze oscillations with different ve- 
locities and frequencies, which make amplitudes in the 
range of several tens of nanometers. We assume that 
the expectation value of the qubit z— coordinate follows 
a saw-like function, 



<*(*)> = 2t 
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where |_- • -J stands for the floor function, T is the pe- 
riod, and £ = vT/2 is the maximum displacement 
from the initial position. Every time the electron re- 
turns to the initial position, once per period, it scans 
exactly the same effective field. As a consequence, 
the spin dynamics experiences strong correlations which 
yield the memory effects. In the oscillatory regime, 
the correlator for the precession rate satisfies the con- 
ditions ((Sli(i)Sli(T(l2t/T\+l)-i))) = and 
((Qi(t)Slilt + nT))) = ((Of)), where n > is an inte- 
ger. 

We begin with the regime where TT V is of the order of 1 
and the typical electron spin strongly changes in a single 
oscillation. The spin polarization behavior is shown in 
Fig. 2J We begin with the period T = 7.5 ns and velocity 
v = 20 nm/ns, where spin relaxation during half a period 
is not very strong yet. As a result, the polarization first 
decays and then increases to a stable value with tiny saw- 
like peaks on the top. For a twice larger period, T = 15 
ns and v = 20 nm/ns, the product TT V is larger by a fac- 
tor of 2, spin relaxation per period is considerably more 
efficient, and the calculated dynamics already presents 
clear periodic peaks above the reached plateau. For a 
yet smaller frequency and speed, T = 30 ns and v = 10 
nm/ns, one observes an almost flat stage followed by a 
series of sharp revivals periodically repeated in time, see 
Fig. [4J Here, in the first period, the polarization relaxes 
down to zero. Throughout the second period, the qubits 
show a revival to 1/3 of the initial spin. Beyond two pe- 
riods, relaxation and memory effects both contribute to 
yield a periodic picture of revivals up to one third and 
decays down to zero. 

In the next case, we consider very high velocity and 
frequency, TT V <C 1, resulting in a small spin precession 
angle in one period. When v = 160 nm/ns and T = 2.5 
ns, the average spin polarization decays to a minimum 
of about 10% of the initial value and monotonically in- 
creases up to a constant plateau of 1/3; the ensemble 
dephasing time is roughly 14 ns, as can be seen in Fig. [5] 
Due to the extremely high speed, the qubit spin cannot 
follow, by precessional motion, the fast changes in the ef- 
fective nuclear magnetic field. The characteristic behav- 
ior is given in Fig.[5j similar to that depicted in Fig. [2] for 
static qubits. Here, the behavior of polarization in time 
cannot be modified by tuning the qubit width. 

In the oscillatory regime, the relaxation part of the 
time-dependence is Gaussian rather than exponential 
32[ 33| , as we explain below. At any t, ensemble mean 
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FIG. 5: (Color online) Average polarization of an oscillating 
qubit when T — 2.5 ns, v — 160 nm/ns, and £ = 10 nm. The 
inset shows details of the short-term behavior changing slope 
every half a period shown by vertical dotted lines. 



value of spin z— component is ((cos# z ))/2, where Z 
is the corresponding precession angle; this mean value 
is ((cos0 z )) = exp(-y/7rT v t) for < t < T/2, and is 
multiplied by factor n, ((cos0 z )} = exp(— y / 7rn 2 r v t), at 
t = nT/2. Thus, in this regime, a new spin relaxation 
rate is established, r mem = y/Y v /T > r„, similar to the 
memory effects in the systems with random s pin- orbit 
coupling subject to external magnetic field (32I. l34j . On 
the other hand, since r mem = r st ^r^jT <C r st , this 
relaxation is found even slower than for static qubits. 
Moreover, due to the memory effects, the slopes of 
]n(((s z (i)))/((s z (0)))) at t = nT/2 should change follow- 
ing the ratios (2n + l)/(2n — 1), which are indeed very 
close to those obtained from our numerical data, see inset 
of Fig. 03 

To summarize, we have investigated the spin dynam- 
ics due to hyperfine coupling of electrons embedded in 
quantum dots propagating along quantum wires. This 
behavior of spins is qualitatively different from the typi- 
cal single-parameter time dependence in static dots. For 
straightforward motion with constant velocity, the spin 
relaxation is close to exponential towards zero and, due 
to motional narrowing, its rate does not depend on the 
spatial width of the packets but solely on their speed. On 
the contrary, in the oscillatory regime, the decay can be 
Gaussian rather than exponential, and the polarization 
revives afterwards up to one third of the initial value. 
Two modes characterize the oscillatory regime: The spin 
revival can be reached through periodic peaks and val- 
leys, and, in the other mode, through a monotonous 
increase of polarization towards a stable plateau. The 
shown dynamic trends, especially interesting for oscil- 
lating qubits, adds to the knowledge of the behavior 
of "flying" electrons in semiconductor nanowires. They 
also broaden the applicability of such nanosystems in the 
fields of spintronics and quantum information technology. 
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